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Abstract 

The effect of Pd doping on the superconducting phase diagram of the unconventional 
superconductor UPt 3 has been measured by (magneto Resistance, specific heat, thermal 
expansion and magnetostriction. Experiments on single- and poly crystalline U(Pti. x Pd x )3 for 
x< 0.006 show that the superconducting transition temperatures of the A phase, T* , and of 
the B phase, T~ both decrease, while the splitting AT C increases at a rate of 0.30+0.02 
K/at.%Pd. The B phase is suppressed first, near x= 0.004, while the A phase survives till 
x= 0.007. We find that AT c (x) correlates with an increase of the weak magnetic moment m(x) 
upon Pd doping. This provides further evidence for Ginzburg-Landau scenarios with 
magnetism as the symmetry breaking field, i.e. the 2D E-representation and the ID odd parity 
model. Only for small splittings AT c °c m 2 ( T c + ) (AT C < 0.05 K) as predicted. The results at 
larger splittings call for Ginzburg-Landau expansions beyond 4th order. The tetracritical point 
in the B-T plane persists till at least x= 0.002 for B_L c, while it is rapidly suppressed for Bll c. 
Upon alloying the A and B phases gain stability at the expense of the C phase. 
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1. Introduction 

The superconducting instability in heavy-electron compounds [1,2] continues to attract a great 
deal of attention. In the past years much research has been directed towards the close 
connection between superconductivity and magnetism in heavy-electron materials [3]. The 
principle research issues which have emerged are: (i) spin-fluctuation versus phonon mediated 
superconductivity, (ii) the symmetry of the superconducting gap function, and (iii) the 
interplay of magnetic order and superconductivity. Among the heavy-fermion superconductors 
UPt3 with a superconducting transition temperature T c ~ 0.55 K [4] is regarded as exemplary. 
The low-temperature normal state is characterized by pronounced antiferromagnetic spin- 
fluctuation phenomena (T*~ 20 K) and incipient magnetism [5], which give rise to a strong 
renormalization of the effective mass, of the order of 100 times the free electron mass. 
Neutron-diffraction experiments have shown that superconductivity in UPt3 coexists with 
antiferromagnetic order, which develops below a Neel temperature T N ~ 6 K [6]. The 
antiferromagnetic order is unconventional in the sense that the ordered moment squared m (T) 
grows quasi-linearly with temperature. Moreover, the size of the ordered moment is extremely 
small m= 0.02±0.01 (iB/U-atom. The superconducting ground state is difficult to reconcile 
with strong magnetic interactions and, therefore, it has been suggested that superconductivity 
is mediated by antiferromagnetic interactions rather than by phonons [7]. However, decisive 
experimental evidence for this is still lacking. More recently, it has been argued that 
superconductivity is a more general property of heavy-fermion antiferromagnets close to a 
quantum critical point [8]. In the case of UPt3 the quantum critical point might be reached by 
doping [5], but the concurrent non-Fermi-liquid behavior has not been signaled so far. In the 
past decade, evidence has accumulated that superconductivity in UPt 3 is truly unconventional, 
i.e. the symmetry of the superconducting gap function is lower than that of the underlying 
Fermi surface [9]. Evidence for this is in part presented by the power- law temperature 
dependence of the electronic excitation spectrum below T c , indicating point nodes and/or line 
nodes in the gap [10]. The discovery of a multicomponent superconducting phase diagram 
with three vortex phases in the field-temperature plane [11-14], and the subsequent 
explanation within the Ginzburg-Landau theory of second order phase transitions (see Ref. 15 
and references therein) is in general considered as hard proof for unconventional 
superconductivity . 
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UPt 3 is the only known superconductor with three different superconducting vortex 
phases. In zero magnetic field two superconducting phases are found, the A phase below 
T c + = 0.54 K and the B phase below T c " = 0.48 K. In a magnetic field the A phase is 

suppressed, while the B phase transforms into a third phase, labeled C. The three phases meet 
in a tetracritical point. The phenomenology of the phase diagram has been studied extensively 
using Ginzburg-Landau (GL) theory, where the free energy functional is derived exclusively 
by symmetry arguments (the symmetry group for UPt 3 is D 6h ). A number of GL models have 
been proposed [15-22] in order to explain the zero-field splitting AT C = T c + -T c " [11] and the 

topology of the phase diagram in magnetic field [12-14] or under pressure [23, 24]. Most of 
the GL models require an unconventional superconducting order parameter. The most 
plausible GL models which have been worked out to understand the phase diagram of UPt 3 
fall into three categories: (i) the degeneracy of a two-dimensional (2D) even or odd parity 
order parameter is lifted by a symmetry breaking field [16-19], (ii) the spin degeneracy of a 
one-dimensional (ID) odd parity order parameter is lifted by a symmetry breaking field under 
the assumption of a weak spin-orbit coupling [15, 20], and (iii) there is an accidental 
degeneracy of two nearly degenerate ID representations [19, 21]. However, no consensus has 
been reached, as each of the three models only partially describes the field and pressure 
variation of the superconducting phases. As regards the first two scenarios a key issue is to 
identify the symmetry breaking field (SBF). Experimental evidence that the anomalous weak 
antiferromagnetic order which sets in at T N ~ 6 K acts as the SBF is at hand [25]. Another 
candidate for the SBF is the incommensurate structural modulation which has been detected 
by transmission electron microscopy [26]. However, its precise role remains unexplored. 

In this paper we focus on the GL models with the degeneracy lifted by a SBF [15-20]. 
More specifically we investigate the role of the small-moment magnetism as SBF. Within the 
model (see section 2), the splitting of the superconducting transition temperature is 

proportional to the strength of the SBF or A7>« £, where m . Direct evidence for the 

coupling between AT C and m was deduced from specific-heat [23] and neutron-diffraction 

[25] experiments under hydrostatic pressure. It was observed that both AT Z , determined by 
specific heat, and m (T c ), measured by neutron diffraction under pressure, vary linearly with 
pressure and vanish at a critical pressure p c ~ 3 kbar. We utilize another route to verify the 

coupling between AT C and m 2 , namely by doping UPt3 with small amounts of Pd. 
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Vorenkamp and co-workers carried out specific-heat experiments on poly crystalline 

samples of U(Pti_^Pd^)3 (x< 0.002) and showed that AT C almost doubles with respect to pure 
UPt 3 for the x= 0.002 compound [27]. This then directly prompted the question whether the 

enhancement of AT C is due to the increase of the ordered moment m. Since it was known that 
for 0.02oc<0.07 pronounced phase-transition anomalies in the thermal and transport 
properties signal an antiferromagnetic phase transition of the spin-density-wave type [5], we 
conducted a neutron-diffraction study on single-crystalline samples in order to investigate m 

as function of Pd content over a wider range of x, including the region where AT C is observed 
to increase. These results are reported in Ref. 28 and the conclusions are as follows. The 
small-moment antiferromagnetic order (SMAF) is robust upon doping with Pd and persists 
until at least x= 0.005. The ordered moment grows from 0.018±0.002 (i B /U-atom for pure 
UPt 3 to 0.048±0.008 (l B /U-atom for x= 0.005. For the SMAF T N ~ 6 K and does not vary with 
Pd contents. For x> 0.01 a second antiferromagnetic phase is found, for which at optimum 
doping (x= 0.05) T N attains a maximum value of 5.8 K and the ordered moment equals 
0.63±0.05 (l B /U-atom. For this large-moment antiferromagnetic order (LMAF) T N (x) follows a 
Doniach-type diagram. From this diagram it is inferred that the antiferromagnetic instability in 
U(Pti. x Pd x ) 3 is located in the range 0.5-1.0 at.% Pd. 

In this paper we present a study of the superconducting properties of U(Pti_ x Pd x )3. The 

main objectives of this work are: (i) to determine AT c (x) by means of specific-heat 

experiments, (ii) to test the SBF model by relating AT c (x) to the ordered moment mix), and 
(iii) to investigate the effect of Pd doping on the superconducting phase diagram in the B-T 
plane by means of magnetotransport and dilatometry experiments. The paper is organized as 
follows. In section 2 we review the basic relations for the SBF scenario. In section 3 we 
concentrate on the sample preparation process and the characterization of the samples by 
means of electrical resistivity. In section 4 we present and analyze the specific-heat of 
U(Pti-^Pd^)3 in the vicinity of the double superconducting transition. In sections 5 and 6 we 
present the magnetoresistance, thermal expansion and magnetostriction data and in section 7 
we construct the phase diagrams in the B-T plane for x= 0.002. In section 8, we extract the 
Ginzburg-Landau parameters, while the SBF model is tested in section 9. Finally, we present 
the concluding remarks in section 10. Parts of these results have been reported in a 
preliminary form in Refs. 29-31. 
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2. The SBF scenario 

The SBF scenarios can be divided into two categories: (i) the degeneracy of a 2D even or odd 
parity order parameter is lifted by a SBF [16-19] and (ii) the spin degeneracy of a ID odd 
parity order parameter is lifted by a SBF under the assumption of weak spin-orbit coupling 
[15, 20]. The irreducible representations for the superconducting gap with the appropriate D^h 
symmetry of UPt 3 have been tabulated by e.g. Yip and Garg (Ref. 32). We first concentrate on 
the 2D representation called the E-representation model. For a 2D representation with even 
parity, E\ g or E 2g , or odd parity, E\ u or E 2u , the superconducting gap function is given by 
Ag(k) = T| x r£ jX (k) + T^r^k), where Te,x and T^y are the basis functions for the relevant 2D- 
representation. The complex vector r|= (r| x , r| y )= (lr| x le 1<Px , lr| y le 1<Py ) determines the order 
parameter. The free energy functional can be written as the sum of three terms [15-20]: 

F = Fs + Fm+Fsm (1) 
Here F s is the free energy functional of the superconductor 

F s =a s (r-r c )|Ti| 2 +|(3 1 |Ti| 4 +|(3 2 |Ti 2 | 2 (2) 

where the coefficients as, (3i and P2 are stability parameters. The contribution from the 
magnetic order to the free energy is given by 

F M =a M (r-r N )m 2 +|(3 M m 4 (3) 

where m= (m,0,0) is the small ordered moment oriented along a principal axis in the basal 
plane (T< T N ) and oc M and (3 M are stability parameters. The mixing term of magnetic order and 
superconductivity can be written as [16-19]: 

^ M =-7^ 2 (r| x -r| y ) (4) 
where 8= ym is the symmetry breaking field. By minimizing the free energy it follows that the 
single superconducting transition at T c splits into two transitions at T c + and T~ , where 

AT c =T;-T c -= 8 ^^^m 2 (5) 

Here g= 1 for the E-model. The ratio (Pi+f^Vfb can be determined from the measured step 
sizes in the specific heat at T c + and T~ : 
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Ac(r c -)/r; =1| p 2 
Ac(t c + )/t; p 1 



(6) 



Here the steps in the specific heat are measured relative to the normal state. The weak- 
coupling value for pV|3i is 0.5. 

The ID representation model with odd parity yields very similar expressions. The three 
component order parameter is r|= (r| x , r| y , r] z )= (lr] x le 1(px , lr| y le 1(Py , lr| z le 1<Pz ) and the gap function is 
given by [20] : 



with %x=io y Ox, where the o's denote the Pauli spin matrices. The complex coefficients x\x are 
characterized by a spin index X. The orbital part, /(k), belongs to the ID representation A2 U , 
Bin or B2u- The free energy functional is expressed as in the E-model using equations 1-3. The 
coupling term of the magnetic and the superconducting order parameter consists of three 
components and equation 4 now reads: 



For the ID model AT C is given by equation 5 with g= 3/2, while the expression for pV|3i is the 
same as in the E-model (Eq. 6). Note that in Ref. 20 an incorrect expression is given for T~ in 

which (3i and are interchanged in the numerator. 

As it is our purpose to verify Eq. 5 by experiments, one also needs, besides values for AT C 
and pV|3i, which can be deduced from the specific-heat data, and the value of m, which 
follows from the neutron-diffraction experiments [28], an estimate for the parameter lyl/cxs- 
We comment on this in section 9. 

3. Experimental 

The data reported in this paper have been taken on annealed poly crystalline (x< 0.006) and 
single-crystalline samples (x< 0.002). Poly crystalline material was prepared by arc-melting the 
constituents in a stoichiometric ratio in an arc furnace on a water-cooled copper crucible under 
a continuously Ti-gettered argon atmosphere (0.5 bar). As starting materials we used uranium 
(JRC-EC, Geel) with a purity of 99.98%, and platinum and palladium (Johnson Matthey) with 
purity 5N. Poly crystalline material with low Pd contents (x< 0.01) was prepared by using 
master alloys (e.g. 5 at.% Pd). Single-crystals with x= and 0.002 were prepared in a mirror 



A(k)= £%/(k)x, 



(7) 



X=x,y,z 




(8) 
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furnace (NEC-NSC35) using the horizontal floating zone method. A single-crystalline sample 
with x= 0.001, was pulled from the melt using a modified Czochralski technique in a tri-arc 
furnace under a continuously Ti-gettered argon atmosphere. For annealing, the samples were 
wrapped in tantalum foil and put in a water-free quartz tube together with a piece of uranium 
that served as a getter. After evacuating (p< 10" 6 mbar) and sealing the tube, the samples were 
annealed at 950 °C during 4 days. Next the samples were slowly cooled in 3 days to room 
temperature. Several samples were investigated by Electron Probe Micro Analysis (EPMA), 
but the concentration of Pd is too small to arrive at a quantitative composition analysis. In the 
following, the value of x is the nominal composition. Samples with appropriate dimensions 
and weights were obtained by means of spark erosion. 

In order to characterize the prepared materials the electrical resistivity, p(7), was 
measured on bar-shaped samples. The results for the poly crystalline samples (x= 0, 0.0025, 
0.003, 0.0035, 0.004 and 0.005) are reported in Ref. 30. The data above T c + are well-described 

by the Fermi liquid expression p= po+AT 2 (T< 1 K). The residual resistivity, po, is extracted by 
extrapolating the AT 2 term to T= 0. For x= 0, the residual resistance ratio 
RRR= R(300K)/R(0)~ 1000, indicating a high quality of the pure compound, while 
T c + = 0.57 K. Upon alloying, p increases linearly with x at a rate of 9.6±0.2 (iQcm/at.%Pd, 

which ensures that palladium dissolves homogeneously in the matrix. Also T c + varies 
smoothly with Pd content and the critical concentration for the suppression of 
superconductivity is x c , S c~ 0.007. In case of the single-crystalline materials, p(7) was obtained 
for a current, /, along the a and c-axis. The residual resistivity, po, a , amounts to 0.52, 1.6 and 
2.5|lQcm, while p , c , amounts to 0.18, 0.75 and 1.02 (iQcrn, for x=0, 0.001 and 0.002, 
respectively. For pure UPt 3 we obtain RRR values of -460 and -720 for III a and III c, 
respectively. T c + is suppressed at a rate 0.77 K/at.%Pd. In the following section we compare 

the resistively determined T c + with the bulk value determined by the specific heat. 

4. Specific heat of U(Pt!. x Pd x ) 3 

The specific heat, c(T), of a series of U(Pti_ x Pd x ) 3 samples was measured using the relaxation 
technique. Experiments have been carried out on annealed single-crystalline samples with 
x= 0.000, 0.001 and 0.002 and on annealed poly crystalline samples with x= 0.000, 0.0025, 
0.003, and 0.004. The typical sample mass was 80 mg. The results are shown in Fig. 1 in a 
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plot of c/T versus T. At least three interesting features strike the eye: (i) T c + and T c " are well 
resolved for x< 0.003, while for x= 0.004 only T c + is resolved (T> 0.1 K), (ii) both T c + and T~ 

decrease smoothly with Pd concentration, while Ar c increases, and (iii) the overall height of 
the jumps in c/T at T c + and T~ decreases with increasing x. The results for the single- 
crystalline samples (x< 0.002) are in good agreement with those obtained by Vorenkamp et al. 
[27] on poly crystalline material. In order to determine the ideal values for the jumps in the 
specific heat, we have made use of an equal entropy construction at the NA and AB phase 
boundaries. The ideal transitions are represented by the solid lines in Fig. 1. The resulting 

values of T c + , T~ , AT C , A NA c(T c + ) / T c + , A NB c(T~) / T~ , A AB c(T~ ) / T~ and |3 2 /p\ are 

collected in Table I. Here the subscripts NA and NB refer to the step sizes measured with 
respect to the normal state c/T value, while the subscript AB refers to the step size measured 
with respect to the c/T value in the A phase. Below T' c c s (7)= joT+dT 2 , down to the lowest T 

measured (0.1 K). The dT 2 term shows that the superconducting gap function has a line node 
[9] . For T—> K considerable residual yo values are observed which is attributed to impurity 
broadening of the line node [33]. Just as is the case for pure UPt 3 [11,12], Fig. 1 shows that 
the superconducting state entropy exceeds the entropy of the normal state (assuming cn= JnT). 
The extrapolated entropy unbalance for < T < T* is slightly sample dependent in 

U(Pti_ x Pd x )3 and ranges from 6 to 12% of the normal state entropy. The entropy discrepancy 
can be resolved by either an increase of c^/T or a decrease of cJT below 0.1 K. The most 
plausible explanation for the entropy imbalance is offered by the presence of an anomaly at 
0.018 K in the normal state specific heat [34]. The entropy balance is fulfilled when this peak 
is included. 

In Fig. 2 T c + , T c " and AT C are plotted as function of Pd concentration. Both, T* and T~ 
decrease with increasing Pd concentration, but with different rates: dT* /dx= 
-0.79+0.04 K/at.%Pd and dr c "/cU= -1.08+0.06 K/at.%Pd, and as a result AT C increases at a 

rate dAr c /dx= 0.30+0.02 K/at.%Pd. The value dr c + /cU= -0.79+0.04 K/at.%Pd measured by 
the specific heat is within the experimental error equal to the resistive value -0.77 K/at.%Pd. 

Usually, the ratio pV|3i is calculated from the steps Ac/T at T c + and T c " with respect to the 
normal phase (Eq. 6). However, in order to obtain proper values of P2/P1 one should realize 
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that Eq. 6 is only correct for small values of AT C , which is not the case in the doped samples. 
Therefore we use a slightly different relation for pV|3i given here below. The steps Ac/T are 
derived from the GL free energy by Ac/T= -d 2 F/dT 2 . The thermodynamic step in the specific 
heat at T c " can be written as: 

a ab c(t; )/t;= A NB c(r; ) / t; - a na c(t; ) / t; (9) 

The temperature dependence of cIT at the two phase transitions is shown schematically in Fig. 
3. In the 4th order GL model cIT is constrained to be temperature independent (solid line), 
while the measured behavior shows a clear temperature dependence (dotted line). In fact, 
higher order terms need to be taken into account in order to arrive at a temperature dependent 
c/T. In order to arrive at a proper estimate of pVp\ we use the directly measured step 

A AB c(r c ~) / T~ , instead of Ac/T at T' c with respect to the normal state. This results in: 

p 2 = A AB c(r c -)/r; 
p\ A NA c(r;)/r; 

The values of pVp\, determined from Eq. 10, are listed in Table I, and are close to the weak 
coupling limit 0.5. The ratio pV|3i is within the experimental error independent of Pd 
concentration. Note that in a first analysis of the specific-heat data we used Eq. 6 which led to 
a steady decrease of pV|3i upon Pd doping [29]. 



5. The upper critical field 

In order to investigate the effect of Pd doping on the upper-critical field, B c2 (T), we have 
measured the electrical resistivity in field for single-crystalline U(Pti„ x Pd x )3 with x= 0.001 and 
0.002. These experiments were primarily conducted to investigate the presence of a kink in 
B&iT), which locates the tetracritical point in the multicomponent B-T phase diagram of pure 
UPt 3 . The experiments were carried out on bar-shaped samples with the current along the long 
axis (III a). The samples were cut from the same single-crystalline batch as used for the 
specific-heat (section 4) and neutron-diffraction experiments [28]. B c2 (T) was determined by 
resistivity experiments in a transverse constant magnetic field for Bll c and B_L c (i.e. Bll a*, 
where a is taken at right angles to a and c). In Fig. 4 some typical results are shown for 
x= 0.001 (B_L c). In these low magnetic fields (B< 1.3 T) the magnetoresistance is small (less 
than 1% of po per Tesla). At each applied field T c + was determined by the 50% resistivity 
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criterion, and the width of the superconducting transition, Ar c + , was determined by the 10- 
90% resistivity criterion. The resulting upper-critical field curves for B_L c and Bll c are shown 
in Fig. 5, where both axes have been normalized by dividing by . For comparison we have 

also plotted in Fig. 5 the resistively determined B C 2(T) data of pure UPt3 [35, 36]. 

For B_L c clear kinks in B c2 (T) are observed (Fig. 5b). This strongly suggests that in the Pd 
doped samples (x< 0.002) a tetracritical point is present, as for pure UPt3. Upon doping the 
tetracritical point shifts towards lower temperatures and higher fields, which indicates that the 
A phase becomes more stable. For x= 0.001 T cr = 0.309(8) K and B cr = 0.461(8) T, while for 
x= 0.002 T CI = 0.225(8) K and B C1 = 0.490(8) (B_L c). Thus for Bl c the phase diagrams for 
U(Pti„ x Pd x )3 (x< 0.002) have the same topology. For Bll c no distinct anomalies are observed in 
B C 2(T) of the Pd doped samples (Fig. 5a). However, for pure UPt 3 a weak kink in B c2 (T) was 
reported [36], locating the tetracritical point at T cx = 0.45(2) K and B cr = 0.60(2) T. In the 
following section we study the phase diagrams for B_L c and Bll c in more detail by 
dilatometry. 

6. Thermal expansion and magnetostriction 

In order to determine the superconducting phase diagram of the x= 0.002 compound 
dilatometry experiments (thermal expansion and magnetostriction) have been performed. The 
results will be compared with the dilatation experiments on pure UPt 3 reported by Van Dijk et 
al. [37,38]. 

6.1 Experimental 

The U(Pto.998Pdo.oo2)3 sample used for the thermal expansion and magnetostriction 
experiments was cut from the same single-crystalline batch as prepared for the resistivity, 
specific heat and neutron-diffraction experiments. The approximate dimensions of the sample 
along the a, a* and c-axis are 4.0, 5.0 and 3.4 mm, respectively. The sample was mounted in a 
capacitance dilatation cell machined of oxygen-free high-conductivity copper. Two Ru02 
resistors which served as heater and thermometer were glued onto the sample. Length changes 
along the c-axis of the sample were determined with the three-terminal capacitor technique, 
using an Andeen-Hagerling bridge (model 2500E). The sensitivity of the experimental set-up 

o 

amounts to 0.01 A. The dilatation cell was attached to the cold finger of a dilution refrigerator. 
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The coefficient of linear thermal expansion, oc= L~ dL/dT, was measured using a modulation 
technique (f= 0.003 Hz, Ar=5-10mK). The linear magnetostriction, k= (L(5)-L(0))/L(0), 
was measured by sweeping the magnetic field at a relatively low rate (d5/d£< 0.03 T/min) 
while monitoring the length, L, of the sample. The magnetostriction was measured for a field 
along the dilatation direction (Bll c) and at right angles (Bll a). 

6.2 Thermal expansion 

The zero-field temperature variation of the coefficient of linear thermal expansion along the c- 
axis, a c (7), of U(Pt .998Pdo.oo2)3 is shown in Fig. 6, with as inset a c (J)/T. Just as for pure UPt 3 , 
aJT is constant in the normal state, while two clear steps of opposite sign (most pronounced 
in a c (7)/7) mark the double superconducting transition. The superconducting transition 
temperatures have been determined using an equal-length construction and the idealized 
transition is given by the solid line in Fig. 6. For 0.002, T*= 0.381(2) K and 

T~= 0.271(4) K. These values are in excellent agreement with the transition temperatures 

T c + = 0.384(3) K and T~ = 0.276(4) K determined by the specific heat (see section 4). 

However, the value of T* = 0.420(3) K determined resistively is slightly higher. This has also 

been noticed for pure UPt3 [12]. The resistive transition temperature marks the onset of the 
bulk transitions measured by the specific heat and thermal expansion. The difference between 
the resistive and bulk transition decreases in an applied magnetic field. In Fig. 7, a few 
exemplary a c (7)-curves are shown in a constant magnetic field (Bll c and Bll a). Both T c + and 

T~ are suppressed with field, but T c + is suppressed more rapidly than T~ , so that they merge 

at a critical field, 5 cr . The field dependence is anisotropic. For Bll a the transitions merge in 
the field range 0.5-0.6 T, while for Bll c the transitions do merge at about 1.0 T, which is close 
to B C 2 at our lowest temperature (0.075 K). 

Combining the thermal-expansion and the specific-heat data we can determine the 
uniaxial pressure dependence of the superconducting phase transitions using the Ehrenfest 
relation: 

dT = y m Aoc, 
d Pi A(c/T) 

Here p, is the uniaxial pressure along a specific crystallographic axis (i= a, a*, c) and 

CO 

V m = 4.24x10" m /mol is the molar volume. With help of the thermal-expansion data of Fig. 6 
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and the specific-heat steps listed in Table I, we calculate: dT c + /dp c = -0.14(1) K/GPa and 
dT~/dp c = 0.06(1) K/GPa. Thus for uniaxial pressure along the c-axis the splitting, 

Ar c = T c + - T~ , decreases at a rate dA7Vdp c = -0.20(2) K/GPa. These calculated pressure 
dependencies are similar to those determined from the pressure dependence of the specific 
heat for pure UPt 3 , where dT c + /dp c = -0.13(3) K/GPa, d T~/dp c = 0.09(3) K/GPa and 

dATydpc = -0.22(6) K/GPa [39]. Assuming a linear pressure dependence of T c + and T~ , the 

A phase vanishes at p cr = 0.54 GPa and T cr = 0.308 K for x= 0.002, while for pure UPt3 p cr is 
only 0.25 GPa, because of the much smaller zero-pressure splitting, and T cx = 0.459 K. 

6.3 Magnetostriction 

The linear magnetostriction along the c-axis, X C (B), at T= 0.075 K is shown in Fig. 8 for fields 
up to 2 T (Bll c and Bll a). In addition to the normal state contribution to X C (B), a 
superconducting contribution is present below B C 2. For Bll c the normal-state magnetostriction 
(B< 2 T) is well described by a quadratic field dependence, X C (B)= X c (0)+b c B 2 . The coefficient 

7 2 

of the quadratic term, b c , is slightly temperature dependent and is -1.66x10" T" at 
T= 0.075 K. In Fig. 9 we show X C (B) with Bll a at several temperatures. Here the normal-state 

2 7 2 

magnetostriction also follows a B dependence (B< 2 T) with b c is -0.46x10" T" at 
T= 0.075 K. For Bll a the upper critical field B C 2 is difficult to distinguish, while the B-to-C 
phase transition at 5* is visible as a clear kink in the data. For Bll c the situation is reversed: 

B C 2 shows up as a clear anomaly on the X c (5) curve, while B* c does not. For this field 

direction, the superconducting signal, obtained after subtracting the quadratic background 
contribution, is show in Fig. 10. Although the behavior observed for x= 0.002 is in many 
aspects similar to the behavior for pure UPt 3 , two important differences should be noted: (i) 
for x= 0.002 only B C 2 is resolved from the magnetostriction curves for Bll c, while for pure 
UPt 3 both B c2 and B c are resolved, and (ii) for x= 0.002 a significant hysteresis is observed 

for Bll c, which was absent in the data of pure UPt3. In Fig. 11a typical magnetostriction cycle 
(field sweep up and down) is shown. In Fig. 12 we show for both field orientations the amount 
of hysteresis, obtained after subtracting the sweep-up signal from the sweep-down signal. For 
Bll a the hysteresis is negligible, while for Bll c the resulting curve has two peaks. The peak 
just below B C 2 is reminiscent of the peak effect observed in metallic alloys with strong pinning 
of the flux-line lattice. Recently, the peak effect was found in the magnetization of several 
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pure UPt 3 samples [40,41]. The peak effect is expected to become more pronounced upon 
introducing additional pinning centers, e.g. by doping with Pd. The larger low-field peak, 
observed for x= 0.002 at B~ 0.1 T, which is most pronounced for Bll c (see Fig. 12a), has also 
been reported for pure UPt3. This peak, which has a weak temperature dependence, is not 
directly related to the superconducting properties as it is also present in the normal state. The 
origin of this anomaly remains unclear, but it has been suggested that it is related to a meta- 
stable magnetic state [42] . 

With help of the measured discontinuities at the superconducting transitions the field 
dependence of the transition temperatures can be estimated with the Ehrenfest relation: 



Ax. 

= (12) 

J P Ace,. 



dB 

Here x,= dX,/d5 where i refers to the principal crystallographic directions. The initial field 
dependencies of and T~ are determined by the B= thermal expansion data and by 

extrapolation of Ax c to 5— > 0. The values determined in this way are listed in Table II and 
should be compared to the directly measured slopes of the phase lines of the superconducting 
phase diagram of U(Pto.998Pdo.oo2)3 (see section 8). 



7. Superconducting phase diagram of U(Pt .99sPdo.oo2)3 

By locating the anomalies at the superconducting phase transitions determined by our 
dilatometry experiments in the B-T plane we have constructed the superconducting phase 
diagrams of U(Pto.998Pdo.oo2)3 shown in Fig. 13. The NA transition is detected by both thermal 
expansion and magnetostriction, while the AB transition shows up only in the thermal 
expansion. The NC phase line, which is only found for B_L c, has a very weak signature in the 
thermal-expansion data and was therefore complemented by the B C 2(T) data measured 
resistively (Fig. 5b). In this field range the resistive and bulk T c + are equal within the 

experimental accuracy. For B_L c the tetracritical point is located at T CT = 0.205(4) K and 5 cr = 
0.556(8) T. For Bll c the x= 0.002 compound has no tetracritical point (7> 0.075 K), which 
presents a striking difference with respect to pure UPt3. 

The AB phase line of U(Pt .998Pdo.oo2)3 measured for B_L c shows a clear change of slope 
at B= 0.2 T. For pure UPt3 a similar kink was observed, albeit at a lower field B= 0.1 T. It has 
been suggested that this kink arises from a coupling of the superconducting order parameter to 



13 



Version: 22 March, 1999 



the meta-stable magnetic state [42]. For U(Pt .998Pdo.oo2)3, however, the change of slope does 
not coincide with the low-field anomaly observed in the magnetostriction at B= 0.1 T, and the 
origin remains unclear. 

In Table in we compare the measured slopes of the phase lines near the tetracritical point 
with the calculated ones using the Ehrenfest relation (Eq.12). Within the experimental 
accuracy the data agree, which demonstrates their internal consistency. Near the tetracitical 
point the thermodynamic steps should follow the relation A NA + A AB = A NC + A BC , where A is 

Ac/r, Aoc or Ax. We have checked that this relation holds for Ax. The thermodynamic stability 
of a phase diagram with a tetracritical point, where at least three second order phase-transition 
lines meet, leads to strict conditions for the slopes of the four phase lines as formulated in Ref. 
43. In the case of pure UPt 3 these conditions were satisfied [38]. In order to investigate the 
thermodynamic stability of U(Pto.998Pdo.oo2)3 additional specific-heat measurements in an 
applied magnetic field are needed. 



8. Ginzburg-Landau parameters of U(Pti. x Pd x ) 3 

The temperature derivative of the thermodynamic critical field, dBJdT, near T c , can be 
estimated from the jump in the specific heat at the superconducting transition in zero field 
[44]: 

V fdB V 

A(c/T) = ^- ^ (13) 



Ho 



dT 



Here |io the permeability of free space. The thermodynamic critical field is related to the upper 
critical field by B c2 = 2kB c , where K is the (isotropic) Ginzburg-Landau parameter which 
characterizes the superconducting state [44]. The GL parameter is defined as K= A/£, where X 
is the penetration depth and £ is the coherence length. Since we are dealing with a hexagonal 
strongly anisotropic material K, X and £, are anisotropic. In that case, the upper critical field is 
given by B' c2 = 4> / (271^^) , where $o is the flux quantum and i, j and k are the principal 

crystallographic directions. Assuming that the parameters in the basal plane are isotropic the 
following relations between B C 2 and B c hold: 

c2 V a c c (14) 

fl c c 2 =V2K a fl c 
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Here the superscripts a and c refer to the direction of the magnetic field, K a = X a /^ a and 
Kc= X c /^c- The average GL parameter is defined as K av = (k^k c ) 1/3 . 

We have evaluated the various GL parameters and the temperature derivatives of the 
upper critical fields from the measured data. The results for 5— > are listed in Table IV. The 
value of d5 c /dr for the A phase has been determined using A NA c(r c + )/ r c + . The values of 

d5 C 2/dr determined from the step in the magnetostriction and thermal expansion (see Table II) 
are in reasonable agreement with the values determined directly from the slope of the phase 
diagrams. 

UPt 3 is an extreme type-II superconductor with for pure samples X> 6000 A and %~ 120 A 
so that K~ 50. Upon Pd doping (x< 0.002) K av remains roughly constant, while K a decreases 
and Kc increases (see Table IV). Substituting Pd makes the superconducting properties less 
anisotropic, which is also reflected in the ratio of the anisotropic quasiparticle masses, 
determined by m c l m a = (5 c c 2 / 5 c a 2 ) 2 (see Table IV). The lower-critical field, B cl , is related to 

the thermodynamic critical field according to 5 cl = B c ln(K) / (a/2k) , from which it follows 

that B ci is about 6% of B c . Such small values of B c \ have not been probed in our dilatometry 
experiments. 

9. Testing the SBF model 

One of the main objectives of the specific -heat experiments on the U(Pti_ x Pd x ) 3 system was to 
determine the superconducting splitting AT C as function of x. From the data in Fig. 2 we 
conclude that Ar c increases linearly with x at a rate dAr c /cbc= 0.30+0.02 K/at.%Pd. Within the 
GL models presented in section 2, AT C is proportional to the strength of the SBF or assuming 
that the ordered moment of the SMAF is the SBF AT c <x m 2 (T c + ) (see Eq. 5; we comment on 

the prefactor (l7l/a s )((3i+(32)/(32 at the end of this section). In order to determine m 2 (x) we have 
recently carried out neutron-diffraction experiments [28] on single-crystalline samples. The 
ordered moments at T c + are 0.018(2), 0.024(3), 0.036(3) and 0.048(8) |i B /U-atom, for x= 0, 

0.001, 0.002 and 0.005, respectively. In Fig. 14 we have traced Ar c as a function of m 2 (T c + ) 

for x< 0.002. It is interesting to compare this result with AT C <* m 2 ( T c + ) obtained by the 

hydrostatic pressure experiments [25], because doping increases AT C and hydrostatic pressure 
decreases AT C . A direct comparison is not possible because of the relatively large uncertainty 
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2 2 

in the absolute value of m (p= 0)= 0.03±0.01 (l B /U-atom. Therefore, we scaled the m (p) 
values of Ref. 25 such that m (p= 0)= 0.02 (i B /U-atom. The error bars for the pressure data 
correspond to the relative errors determined by counting statistics, while for the Pd doping 
these are absolute errors. After including the pressure data in Fig. 14, we notice the following 

three points: (i) both the Pd doping and pressure data sets collapse onto one curve, (ii) Ar c is a 

smooth function of m 2 ( T c + ), and (iii) AT C <* m 2 ( ) but in a limited range AT C < 0.05 K. The 

latter result shows that the simple GL models presented in section 2 break down for splittings 

AT C > 0.05 K. This is not unrealistic because the applied Ginzburg-Landau expansion is only 

valid for AT C /T C « 1. Clearly, for enhanced splittings a more sophisticated Ginzburg-Landau 
expansion with terms beyond 4th order is needed. We conclude that there is a clear correlation 

between AT C and m 2 ( T c + ), which is in line with the SMAF acting as SBF. 
10. Concluding remarks 

The experimental results reported in sections 3-6 show that the unconventional 
superconducting properties of UPt3 are extremely sensitive to Pd doping. First of all, 
resistivity experiments show that the A phase signaled by T c + is completely suppressed at a 
critical concentration x c>sc is -0.007. Secondly, the specific-heat experiments show that the B 
phase, marked by the second transition at T c " , is suppressed even more rapidly, with 

x c ~ 0.004. Thirdly, AT C increases with Pd contents. One of the main objectives of the present 
work was to investigate the GL model formulated in section 2, especially with respect to 

whether the SMAF acts as the SBF. Indeed, we find a close correlation between AT c (x) and 

m 2 (T c + ). However, only for AT C < 0.050 K the proportionality between AT C and m 2 , predicted 

by the SBF model, holds. For AT C > 0.050 K m grows more rapidly. The failure of the model 
for larger splittings is attributed to the limited applicability of the simple GL E-representation 

and ID odd parity models. The 4th order expansion near T c is only valid for AT C /T C « 1. 

While SMAF and superconductivity coexist, evidence is accumulating that LMAF and 
superconductivity compete. Recent neutron-diffraction [28] and (iSR [45] experiments 
indicate that the critical concentration for the onset of LMAF is near x c ^= jc CjSC ~ 0.007 [30]. In 
order to put this on firm footing additional |iSR experiments are in progress. The competition 
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between superconductivity and LMAF lends further support for spin-fluctuation mediated 
superconductivity. 

The effects of Pd doping (this work) and hydrostatic pressure [23] on the stability of the 
A phase are opposite. It is interesting to note that this also holds for the B and C phases 
[24, 37, 38], which is most clearly observed for B_L c. By applying hydrostatic pressure the 
tetracritical point in the B-T plane shifts to lower fields. Upon increasing pressure first the A 
phase disappears (at p c ~ 0.35 GPa), followed by the B phase, so that the C phase is the most 
stable phase under pressure [24, 37, 38]. For Pd doping the contrary takes place. Upon doping 
the tetracritical point shifts to higher fields, and the A phase gains stability at the expense of 
the B and C phases. Note that the C phase is completely suppressed for U(Pto.99sPdo.oo2)3 in the 
case Bll c. The normal-state properties of UPt 3 react upon Pd doping also in an opposite way 
to hydrostatic pressure. Experiments demonstrate that doping of 1 at.% Pd corresponds to an 
external pressure of about -0.33 GPa [46, 47]. This illustrates that the change of the normal- 
state properties is not governed by the volume, because both Pd doping and applying pressure 
reduce the unit cell volume. Instead, these changes can be explained, to a certain extent, by the 

change in the c/a ratio. In the case of Pd substitution A(c/a)/(c/a)= -0.6xl0" 4 per at.% Pd, while 
for hydrostatic pressure, because of the anisotropic compressibility (k c < K a ), 

A(c/a)/(c/a)= 1.3xl0" 4 per GPa [5], hence, doping 1 at.%Pd corresponds to an applied 
hydrostatic pressure of -0.5(2) GPa. For the stability range of the A phase we do not arrive at 

the same numbers. For pure UPt 3 , dA7Vdp= -0.19 K/GPa [23], while dA7ydx= 
0.30 K/at.%Pd. Thus in this case 1 at.%Pd corresponds to -1.6 GPa. 

In analyzing the specific-heat data around the double superconducting transition, we have 
provided evidence that magnetism provides the SBF. This restricts the choice of the GL 
models to the E-representation model, which applies for both even and odd parity states, and 
to the ID odd parity model. The latter model relies on a weak spin-orbit coupling. In zero 
magnetic field both models give identical results, but they differ in predicting the field and 
pressure dependence of the superconducting phases. Notably, a tetracritical point for all field 
directions is only possible in the E-model under certain conditions and certain symmetries 
(Ei g [48] or E 2u [15]), while no additional constraints are needed in the odd-parity ID model. 
As regards, the pressure dependence, the E model predicts the B phase to be the stable phase 
under pressure. A recent refinement of the odd-parity ID model shows that the C phase is 
most stable under pressure [49]. This is in line with recent pressure studies [24] and 
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dilatometry experiments [38]. Moreover, NMR experiments [50,51] have demonstrated 
convincingly that (i) the Knight shift does not change through the normal-superconducting 
phase transition, and (ii) the effective spin orbit coupling is weak. All these studies provide a 
strong case for the odd parity ID GL model. It is interesting to note that in the refined ID odd 
parity model [49] the antiferromagnetic moment is not static but fluctuates in time. This is 
consistent with recent NMR [50], (iSR [45,52] and neutron diffraction [28] experiments. 

In summary we have studied the superconducting phase diagram of U(Pti_ x Pd x ) 3 , by 
(magneto)resistance, specific heat and dilatometry. Our results in zero field show a strong 

increase of the splitting AT C as function of Pd concentration. AT c (x) correlates with an increase 
of the magnetic moment m(x) upon Pd doping. This provides further evidence for the 
Ginzburg-Landau scenario with magnetism as the symmetry breaking field. The tetracritical 
point in the B-T plane is robust upon alloying for B_L c, at least till x= 0.002, while it is rapidly 
suppressed for Bll c. In a magnetic field the A and B phases gain stability at the expense of the 
C phase upon alloying. In this sense Pd doping and the effect of an external pressure are 
complementary. 
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Table I Parameters deduced from the specific heat of U(Pti_ x Pd x )3. The ratio P2/P1 is 
calculated with help of Eq. 10. The superscripts s and p refer to single- and 



poly crystalline samples respectively. 





t; 


AT C 


A NA c(r c + )/r; 


A NB c(r c -)/r; 


A AB c(r c -)/r; 


P2/P1 


(%) (K) 


(K) 


(K) 


(J/mol K 2 ) 


(J/mol K 2 ) 


(J/mol K 2 ) 




0.00 p 0.560(3) 


0.506(3) 


0.054(4) 


0.23(1) 


0.34(1) 


0.14(1) 


0.60(6) 


0.00 s 0.543(3) 


0.489(3) 


0.054(4) 


0.26(1) 


0.35(1) 


0.13(1) 


0.50(5) 


0.10 s 0.437(3) 


0.355(3) 


0.082(4) 


0.21(1) 


0.26(1) 


0.12(1) 


0.57(7) 


0.20 s 0.384(3) 


0.276(4) 


0.108(5) 


0.17(1) 


0.19(1) 


0.10(1) 


0.58(8) 


0.25 p 0.362(3) 


0.236(4) 


0.126(5) 


0.18(1) 


0.18(1) 


0.08(1) 


0.44(9) 


0.30 p 0.313(4) 


0.163(5) 


0.150(6) 


0.13(1) 


0.09(1) 


0.06(1) 


0.46(11) 


0.40 p 0.222(5) 






0.07(1) 









Table II Thermodynamic quantities Aoc c and Ax c for x= 0.002 
at the superconducting transitions T c + and T~ in zero 



field. The step Ax c is determined in the limit 5— > 0. 







7T(K) 


Aa c (10" 6 K" 1 ) 


-0.56(2) (NA) 


0.14(1) (AB) 


Ax c (lO'V) 


-0.11(1) (Bll a) 
-0.07(1) (Bll c) 


0.01(1) (Bll a) 
0.01(1) (Bll c) 


-Ax c /Aa c (K/T) 


-0.20(2) (Bll a) 
-0.13(1) (Bll c) 


0.07(7) (Bll a) 
0.07(7) (Bll c) 
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Table HI Thermodynamic quantities for x= 0.002 in the vicinity of the 
tetracritical point (Bll a) at T cr = 0.205(4) K and B cr = 0.556(8) T. 
The phase lines between the A, B, C and N phases are indicated by 
NA, NC, AB and BC. 



d77dfl (K/T) 



Aa c (10" 6 K" 1 ) 



At c (10- 6 T 1 ) 



-Ax c /Aa c (K/T) 



NA 



-0.472(8) 



-0.18(2) 



-0.11(1) 



-0.61(12) 



NC 



-0.217(4) 



-0.02(2) 



0.00(1) 



0.0(5) 



AB 



-0.086(1) 



0.06(1) 



0.00(1) 



0.0(2) 



Determined by the relation A NA + A AB = A NC + A BC 



BC 



1.37(6) 



-0.10 



-0.11(1) 



1.1(1) 



Table IV Slopes of the critical fields for the A phase of L^Pti^Pd^, and the calculated GL 



parameters and effective mass ratio. 



X 


dfl c /dr 


d77dfl c a 2 


d77dfl c c 2 


K a 


Kc 


K a v 


m c /m a 


(%) 


(T/K) 


(K/T) 


(K/T) 










0.0 


-0.087(6) 


-0.241(8)* 


-0.093(4)* 


87(6) 


13(1) 


46(3) 


6.7(6) 


0.1 


-0.078(6) 


-0.250** 


-0.124** 


73(6) 


18(1) 


46(3) 


4.1 


0.2 


-0.071(7) 


-0.258(6) 


-0.155(6) 


64(6) 


23(2) 


46(3) 


2.7(2) 



Data taken from ref. 38. 

Average value of the entries for the x= 0.000 and x= 0.002 
compounds. 



23 



Version: 22 March, 1999 



Figure caption 

Fig. 1 Specific heat divided by T versus T of U(Pti^Pd^) 3 for x= 0.000, 0.001 and 0.002 
(single crystals) and for x= 0.0025, 0.003, and 0.004 (poly crystalline samples). The 
solid lines represent ideal transitions determined from an equal entropy construction. 

Fig. 2 T c + , T c and Ar c of U(Pti. JC Pd ;c )3 as a function of Pd concentration, determined from 

the specific-heat data. The solid and open symbols represent single- and poly- 
crystalline data, respectively. 

Fig. 3 Schematic temperature dependence of cIT at the double superconducting transition. 
The solid line represents cIT calculated from the GL free energy, while the dotted line 
reflects the observed behavior. 

Fig. 4 Resistivity of U(Pto.999Pdo.ooi)3 (HI a) in constant magnetic fields Bll a*, ranging from 
to 1.3 T in steps of 0.1 T. For the most right curve B- T and for the most left 
curve B= 1.3 T. 

Fig. 5 The upper critical field of U(Pti_ x Pd x )3, determined by resistivity (III a), in a plot of 
B c2 IT c as a function of TIT C for (a) Bll c and (b) Bll a*. T c = T c + is 0.547(5) K, 

0.466(5) K and 0.420(5) K for x= 0.000, 0.001 and 0.002, respectively. The arrows 
mark the tetracritical points. The data of pure UPt 3 are taken from Refs. 35 and 36. 

Fig. 6 The coefficient of linear thermal expansion along the c-axis (a c ) of U(Pt .998Pdo.oo2)3- 
The solid line represents the ideal transition determined from an equal-length 
construction. In the inset the data is plotted as a c /T vs. T. The transition temperatures 
are T c + = 0.381(2) and T c ~= 0.271(4). 

Fig. 7 The coefficient of linear thermal expansion along the c-axis (a c ) of U(Pto.99sPdo.oo2)3 
in magnetic fields ranging from to 1 T as indicated, with (a) Bll c and (b) Bll a. The 
curves in field are shifted upwards along the vertical axis for the sake of clarity. 
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Fig. 8 Linear magnetostriction along the c-axis (k c ) of U(Pt .998Pd .oo2)3 at T= 0.075 K for 
Bll a and Bll c. The dotted line for Bll c represents the extrapolated normal state 
magnetostriction (see text). 

Fig. 9 Linear magnetostriction along the c-axis (A, c ) with increasing magnetic field of 
U(Pto.998Pdo.oo2)3 for B|| a. The temperature ranges from 0.10 to 0.35 K as indicated. 
The arrows mark the BC and CN transitions (see text). The curves for T> 0.15 K are 
shifted upwards for the sake of clarity. 

Fig. 10 Linear magnetostriction along the c-axis (k c ) of U(Pto.998Pdo.oo2)3 for increasing Bll c 
at temperatures as indicated. The normal state contribution is subtracted (see text). 

Fig. 11 Linear magnetostriction along the c-axis (k c ) at T= 0.075 K of U(Pt .998Pd .oo2)3 for 
Bll c. The normal state contribution is subtracted. The arrows indicate the sweeps up 
and down of the magnetic field. The inset shows a close-up of the irreversible 
magnetostriction peak just below B C 2. 

Fig. 12 The hysteresis in X c for field sweeps up and down (X c , U p-X c ,down) of U(Pto.99sPdo.oo2)3 
with (a) Bll a and (b) Bll c. The peak effect is observed just below B c2 for Bll c. The 
temperatures range between 0.075 and 0.4 K as indicated. 

Fig. 13 The superconducting phase diagram of U(Pto.998Pdo.oo2)3 for B_L c and Bll c, 
constructed from phase transitions detected in the thermal expansion (O) and 
magnetostriction (•). For B_L c the CN phase transition is determined resistively (A). 

Fig. 14 The variation of the splitting AT C as a function of m 2 ( T c + ) for U(Pti. JC Pd ;c )3 (•) and for 

UPt 3 under pressure (□) [25]. For AT C < 0.05 K A7> m as predicted by the GL 
model (see text). The solid line is to guide the eye. 
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